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The  Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles  was  established  in  1962,  at  the  Ohio 
Agricultural  Experiment  Station,  Wooster,  by  the  Agricultural  Research  Service  of  the  United  States 
Department  of  Agriculture,  in  recognition  of  the  need  for  continued  fundamental  research  on  fine  particle 
behavior  and,  in  particular,  the  criticahty  of  such  studies  to  agriculture. 

In  agriculture,  pesticides  in  both  liquid  droplet  and  sohd  particle  form  are  used  extensively,  but  their 
use  is  frequently  hampered  by  lack  of  optimum  efficiency  and  precision  of  apphcation,  despite  the  best 
efforts  of  manufacturers  and  apphcators.  This  results  in  an  excessive  economic  burden  on  agriculture  in 
the  maintenance  of  the  best  of  quahty  in  agricultural  produce,  as  well  as  an  increased  chance  of  contami- 
nating surrounding  areas.  Problems  of  air  pollution  are  also  of  prime  pertinence  and  concern.  It  is 
desired  that  the  results  of  these  investigations  in  fine  particle  behavior  and  its  inherent  subject  matter 
will  serve  agriculture  and  its  allied  industries,  and  other  industries  with  similar  vexations,  in  the  allevia- 
tion of  these  problems. 
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THEORETICAL  ASPECTS  OF  DROPLET  IMPACT  MEASUREMENTS 


Ross  D.  Brazee  ^ 


The  damped  spring-mass  transdueer  is  considered  from  the  standpoint  of  droplet  impaction 
measurements.  The  impulse  concept  of  droplet  impaction  is  mentioned.  Responses  of  the 
transducer  to  step-function  and  random  inputs  are  analyzed  theoretically. 


INTRODUCTION 


The  hnear  spring-mass  system,  in  the  form  of  a cantilever  beam,  has  found  apphcation  to  the  study  of 
liquid  droplet  impaction.^  In  this  paper,  we  will  consider  some  theoretical  aspects  of  the  measurement 
technique.  It  will  not  affect  the  discussion  to  assume  that  the  mass  of  the  system  remains  constant,  i.e., 
no  hquid  is  accumulated  by  the  impaction  surface. 


DROPLET  MOMENTUM 


The  major  force  causing  the  deflection  of  the  spring-mass  system  will  arise  from  the  summation  of 
the  individual  droplet  impulses.  The  mass  of  liquid  in  each  droplet  will  contribute  an  impulse  dependent 
upon  its  momentum  change  in  the  direction  of  the  transducer  deflection.  This  follows  from  the  definition 
of  the  impulse,  or  time  integral  of  a force  F, 


dt—m 


= U — Uo 


(1) 


where  m is  the  mass,  r is  the  position  vector,  and  Uo  and  u are  the  momentum  vectors  at  t — to  and  t = t, 
respectively. 


We  would  expect  that  the  energy  involved  in  the  momentum  change  would  be  spent  in  shattering  of  the 
droplet,  viscous  dissipation  in  the  fluid  and  surrounding  air,  and  transfer  of  momentum  to  directions  other 
than  the  hne  of  impact.  Also,  turbulence  should  result  in  the  residual  fluid  and  heat  should  appear  through 
various  processes.  A target  such  as  soil  would  take  up  energy  through  deformation  and  breakup. 


RESPONSE  OF  THE  MEASUREMENT  TRANSDUCER  TO  A STEP  FUNCTION  INPUT 

The  behavior  of  the  damped  spring-mass  transducer  may  be  described  by  the  differential  equation 

d'^x  , dx  , , , 

+ = <2) 

where  m and  k are  the  apparent  mass  and  spring  constant,  respectively,  c is  the  damping  coefficient,  and 
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F{t)  is  the  impressed  force.  We  apply  a force  F{t)  in  the  form  of  a unit  step  function,  i.e., 


F(^)=l(^-^o). 

(3) 

It  is  convenient  to  define,  in  the  usual 

manner,  the  undamped  natural  angular  frequency. 

co„=  (klm)'^; 

(4) 

the  critical  damping  coefficient. 

and  the  damping  factor 

Cq  — 2m(On", 

(5) 

C = cicc. 

(6) 

Upon  obtaining  the  general  solution  for  (2)  under  these  definitions  and  the  initial  conditions 

x{0-\-)  = xq  = 0; 

(7) 

i(0  + ) = io  — 0, 

(8) 

it  is  possible  to  distinguish  among  the  cases  ^ < 1,  ^=1,  and  ^ >1. 


Case  I.  ^ < 1 (Light  Damping).  For  ^ < 1,  the  general  solution  becomes 


x{t) 


m(ol 


1 — e ^0^ 


cos(o>n(l  — ^2)^  (t  — ^(J))  + sin  — — I 


(9) 


Hence,  the  displacement  x{t)  of  the  transducer  in  response  to  a step  function  input  may  be  represented 
by  a damped  oscillation  superimposed  upon  a step  function. 

Case  II.  i—1  (Critical  Damping).  If  ^=1,  the  apphcation  of  a hmiting  process  and  L’Hospital’s  rule 
to  the  general  solution  yields  the  result 


l(^  — ?o)  r 'I 

x{t)  = — — |l—  [l  + a)„(^— ^o) 


(10) 


This  case  represents  the  transition  from  the  oscillatory  to  the  nonosciUatory  condition  of  the  transducer. 
An  aperiodic  displacement  is  superimposed  upon  the  step  function. 


Case  III.  ^ > 1 (Heavy  Damping).  In  this  case,  the  general  solution  takes  the  form 

*(0=^^r^(l-e-‘“»''-'o'icosh  ( o.„ ( 1 f ( f - to) ) + , . J ^ sinh  ( 1 f (t  - f o) ) 1 1 
main  I L (t  A)  JJ 


(11) 


An  aperiodic  curve  is  superimposed  upon  the  step  function  displacement. 

Obviously,  the  hghtly  damped  transducer  may  indicate  oscillations  which  are  not  representative  of 
the  input.  This  would  be  particularly  true  for  a sudden  change  in  the  input,  such  as  indicated  by  the  above 
results  for  the  step-function  input.  Accordingly,  it  would  appear  desirable  to  damp  the  transducer  in  the 
critical  to  heavy  range. 

The  transducer  wiU  act  as  a band-pass  filter  with  characteristics  dependent  upon  the  natural  frequency 
and  damping.  Therefore,  it  is  now  pertinent  to  describe  its  response  to  a random  input. 
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RESPONSE  OF  THE  TRANSDUCER  TO  A RANDOM  FORCE 


If  the  transducer  is  subjected  to  a random  force  U(t),  the  displacement  becomes  a random  process 
governed  by  the  stochastic  differential  equation 


X (t)  + 2ojnCX{t)  + o)pr(0  = U(t) . 


(12) 


The  symbols  X{t)  and  X{t)  denote  the  random  acceleration  and  velocity,  respectively.  The  processes 
X{t)  and  U{t)  are  assumed  stationary  and  allowed  to  be  possibly  complex  with  the  expectations 

E{X{t)}  = 0;E{U{t)}  = 0,  (13) 


for  convenience,  and 


E{X{t)X*{t)}=(Ta:^;  E{U{t)UHt)}=CTu\ 


(14) 


Since  X{t)  and  U{t)  are  stationary,  they  may  be  analyzed  into  orthogonal  components  in  the  manner 
suggested  by  Bartlett  {p.  168)^  or  Batchelor  (p.  55) The  generalized  harmonic  analyses  of  X{t)  and 
U{t)  in  terms  of  the  orthogonal  processes  Z(o>)  and  JV((o)  will  be 


[? 

J -00 


X(t)=  e*'-dZ(co), 


(15) 


and 


£/(0=  [?^‘^dJE((o). 

J— 00 

Then  (12)  becomes 

r 00  r 00 

I — 2ia)n^(o -h  (o^)dZ(a))  = I e^^^dJF((o), 

J — 00  J —00 

which  upon  combination  with  its  complex  conjugate  form  gives 

l°°[[(‘>>i-o^r+(2to„o^Q,^]dZ(w)dz*(w)=  f ” [ir((u)rfr*(<o). 

J —<x  J ^ ao  J— 00  J— 00 


(16) 


(17) 


(18) 


Taking  the  expectation  of  both  sides  of  (18)  and  then  interchanging  the  operations  of  averaging  and  inte- 
gration, we  obtain 


r 00  r 00 

[(w2- 0)2)2 + (2conO"w2]d£’{Z(a))Z*(a>)}=  \dE{W{ui)W^{oi)]. 

J —00  J —00 


(19) 


It  is  now  helpful  to  define 


and 


Y{oj)  = E{\Z{oj)\^} 
VM=E{\WM\^} 


(20) 


(21) 


^ M.  S.  Bartlett.  An  Introduction  to  Stochastic  Processes.  312  pp.  London,  1956. 

■•G.  K.  Batchelor.  The  Theory  of  Homogeneous  Turbulence.  (Cambridge  University  Press,  London,  1953.) 
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where  Yi(o)  and  V((o)  are  never-decreasing  functions  such  that 

E{X{t)X*{t)]  = E{\X(t)\^}  = cr:,^=  fdE{Z{w)Z*(a)}=  [dy(tu)  = y(=») -y( -=c), 

J—oo  J—x 


(22) 


and  similarly 

with 

and 


E{\U{t)  |2}  = o-t/2=  dVM  = V(oo)  - F(^) , 


Y(-oo)=0;V(-<^)=0, 


F(oo ) = cr/;  F(oo)  = CTu^- 


(23) 

(24) 

(25) 


We  now  make  the  reasonable  assumption  that  X(t)  and  U(t)  have  continuous  spectra,  with  spectral  functions 
Fx(co)  and  Fuict)}  and  spectral  densities  fx((o)  and/jXcu)for  which 


dFxM  =fx((o)da); 

dF  t/(<i))  =fu{(i>)  dco. 


Then  it  follows  that 


E{X*{t)X{t-\-T)}  = Wx{T)  =(Tx^px{r)  =(Tx^  fe^^‘^dFx((o)  =(Tx^  fe^^^fx(co)d(o=  fe*™t/y(o>), 

— nn  •'  — nr  •'  — nr 


(26) 

(27) 


(28) 


where  Wx(t)  and  px(r)  are  the  autocovariance  and  autocorrelation  functions,  respectively.  A siniilar 
procedure  for  U(t)  yields 


F{U*(t)U(t  + T)}  = wu(r)  =ct[/^Pu(t)  =(T[;^  [e^^‘Yu(co)dco=  fe^^^dV(co). 

J— X J— X 

Therefore,  on  the  bases  of  (28)  and  (29), 

dY{o))  =<JxYx{(o)d(o; 
dV{(ii)  =<Ji^fu{(ii)d(ji. 

Returning  to  (19),  we  now  find  it  possible  to  write 

[ [ co^)^+  (2(On0^(^^]^M  = f dV(co), 

J — X J — X 

and  consequently, 

f / \ CTi^fu{(0) 

' cr^^[(a.i-co2)2+(2<o„aW]  ’ 

which  is  the  relationship  between  the  input  and  output  spectral  densities. 


(29) 

(30) 

(31) 

(32) 

(33) 


It  is  of  interest  to  consider  also  the  effect  of  a general  hnear  transducer  or  operator  upon  the  input. 
Here  we  simply  foUow  the  method  outhned  by  Bartlett  {pp.  174—175)  ^ by  defining  the  hnear  ojjerator 


X{t)  = L{U{t)}  = ( g{t  — u)U{u)du. 

J —X 


(34) 


^ See  footnote  3. 
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(35) 


In  the  theory  of  generahzed  harmonic  analysis  (Bartlett,  p.  169),  it  is  shown  that  (34)  is  equivalent  to 

X{t)=  e^^^h{a))dWM, 

J— 00 


where  W((o)  is  an  orthogonal  process  as  in  (16),  and 


For  the  autocovariance  of  X(t)  we  have 


E{X*{t)X{t  + r)}  = 


e'™^(co)/i*(a))c?F(a>), 


(36) 


(37) 


so  that  by  virtue  of  (31), 


o'|/x(a>)  = /^(a>)/^*(cu)o■f/^Xca). 

A particular  case  of  the  linear  operator  L{.  . . .}  is  the  hnear  differential  operator  0{dldt)  =0{D). 
(38)  gives  formally,  then,  for  a 0(i(o)  which  is  non-vanishing  for  real  o> 

^ |0(iO))|2  ’ 

corresponding  to  the  inverse  hnear  operator  L~^.  For  the  damped  oscillator  transducer, 

0(Z))=Z)2-f  \Z>  + 0)2, 


(38) 
Eq. 

(39) 


(40) 


so  that  the  result  given  by  (33)  follows  immediately.  Hence,  it  is  seen  that  the  operational  formalism 
facihtates  the  consideration  of  any  hnear  transducer  for  its  output  as  a result  of  a given  input. 


If  the  spectral  characteristics  of  the  output  are  determined,  then  one  should  be  able  to  estimate  the 
unknown  input  spectrum  by  means  of  (33).  Of  course,  an  alternate  procedure  would  be  to  make  autocor- 
relation measurements  as  indicated  by  (28),  since  the  spectral  and  autocorrelation  functions  are  Fourier- 
transform  related. 


From  (33),  it  is  noted  that  for  light  damping,  the  damped  osciUator  transducer  wiU  have  a response 
sharply  tuned  to  its  natural  frequency.  As  mentioned  in  the  preceding  section,  in  this  state  the  transducer 
wiU  act  as  a band-pass  filter.  It  might  in  this  sense  serve  as  a part  of  a spectrum  analysis  system,  but  not 
necessarily  a convenient  one.  It  would  further,  give  the  false  indications  of  osciUations  on  ordinary 
recorder  traces.  A good  degree  of  damping  would  certainly  be  necessary  for  any  averaging  force  measure- 
ments. Accordingly,  the  value  of  this  particular  transducer  for  more  detailed  studies  of  the  droplet 
impaction  process  may  be  hmited,  or  at  least  subject  to  question. 


IMPLICATIONS  FOR  THE  IMPACTION  PROCESS 


One  may  speculate  about  the  physical  factors  which  would  be  reflected  in  droplet  impaction  spectra. 
It  would  seem  reasonable  to  expect  that  among  these  factors  would  be  found: 

(1)  the  manner  of  breakup  and  dissipation  of  individual  droplet  momenta; 

(2)  droplet  size  spectra; 

(3)  dehvery  characteristics  of  the  droplet  generation  system; 

(4)  characteristic  turbulence  of  the  ambient  medium; 

(5)  effects  of  ambient  medium  turbulence  induced  by  the  droplet  delivery  process. 

(6)  the  nature  of  the  impaction  surface. 


7 


It  would  further  be  expected  that  each  of  these  contributing  factors  would  be  influenced  by  a great  number 
of  elemental  causes.  Hence,  it  appears  that  a most  promising  means  of  an  adequate  description  of  the 
impaction  process  may,  at  least  for  the  present,  be  via  a random  theory  approach.  The  preceding  state- 
ments tacitly  lead  up  to  the  possibihty  that  our  theoretical  considerations  herein  may  be  extendable  beyond 
the  transducer  analysis  to  “real  impaction  processes.”  This  is  not  to  imply  that  a linear  model  wiU  neces- 
sarily be  the  correct  answer  but  perhaps  only  a start.  Stochastic  process  concepts  may  offer  possibihties 
of  forming  generahzations  about  droplet  impaction  processes  and  energy  transfer  which  might  not  otherwise 
be  apparent.  Consequently,  this  mode  of  approach  is  suggested  as  a possible  hne  of  future  research. 
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